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Fidelity approach to quantum phase transitions uses the overlap between ground states of the 
system to gain some information about its quantum phases. Such an overlap is called fidelity. We 
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I. INTRODUCTION 

Quantum phase transitions h^pen when small variations of an external parameter fundamentally change the ground 
state properties of the system [l|, l3- They typically appear when there are competing interactions trying to order the 
sample in different ways and the balance between them is controlled by an external field. 

The fidelity approach to quantum phase transitions uses the overlap between ground states to gain some information 
about the quantum phases [3,0|- For example, we can define fidelity as 

F{9,S) = \{9-S\9 + ^)\, ( 1 ) 

where we assumed that the Hamiltonian H of the system depends on some parameter g (e.g. an external field whose 
variation induces a quantum phase transition), j^) stands for the non-degenerate ground state of the Hamiltonian 
H{9), and ^ is a parameter shift. Besides depending on g and 5, fidelity also depends on TV, i.e., the size of the system 
(number of spins, atoms, etc.). 

The fidelity approach to quantum phase transitions is based on the expectation that fidelity should exhibit a 
marked drop near the critical point, where the ground states of the Hamiltonian are most sensitive to the external 
field. Fidelity has been studied in the following limits: 

• (5 —>■ 0 at the fixed system size N-. one employs here the Taylor expansion to get 

F{g,5) = l-2x{g)5^PO{5^), (2) 

where the linear in 5 term disappears due to the normalization condition {g\g) = 1. Alternatively, one may note 
that such a term could make F{g, <5) > 1 or that it would break the 5 —>■ —<5 symmetry of F{g, 5). The central 
object of this expansion is the fidelity susceptibility xio)- Using the scaling theory of quantum phase transitions 
[I|, Q , it has been shown that at the critical point gc vn & d dimensional system [y, Q 


X(5c) 


( 3 ) 


while far away from it 


x{g) 


N 


\9-9c\ 


2-di' ’ 


( 4 ) 


where v is the universal critical exponent describing the divergence of the correlation length ^(g) near the critical 
point, ^{g) ''^Ig — gd ■ Fidelity susceptibility can be viewed in two ways. First, it can be treated as a function 
whose knowledge leads to the determination of fidelity (HJ through Eq. ©■ One should remember, however, 
that it can be used in such a way only when the term is negligible, which sets conditions on the system 

size N, the distance from the critical point \g — gd, and the field shift 5 . For example, the expansion ([5]) always 
breaks for large enough systems. Second, one can formally define fidelity susceptibility through the relation 


xig) 



dg^ 


( 5 ) 
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and use it to study the quantum phase transitions for any system size N and any parameter g (see e.g. Refs. 
[101 , where hdelity susceptibility gives the diagonal elements of the geometric tensor providing some informa¬ 
tion about the quantum critical points). Finally, we mention that it has been recently proposed that fidelity 
susceptibility can be extracted out from the spectral function [ll| . which can be measured in standard condensed 
matter setups. 

• —>■ oo at the fixed field shift S: this is the limit, where the Anderson orthogonality catastrophe El , i.e., the 

disappearance of the overlap between ground states of the thermodynamically large system, happens. We have 
argued that at the critical point [T^ 


while far away from it 


\nF{gc,S) 

N 


-\ 5 r, 


lnf(g,(5) 


( 6 ) 

( 7 ) 


as long as dv < 2. Eq. ([7]) becomes equivalent to Eq. ([1]) when /\g — gc^ <C 1. We mention that 
— \nF/N is called fidelity per lattice site Ei-El- 

The outline of these lecture notes is the following. Sec. |TT] presents basic concepts associated with the quantum 
Ising model. The main stress in this section is placed on the discussion of the parity of the ground state of the 
Ising chain. In Sec. uni we discuss the exact closed-form expression for fidelity susceptibility of the Ising model. 
Sec. IIVI is devoted to the studies of fidelity in thermodynamically large systems. We show there how the Anderson 
orthogonality catastrophe happens in the Ising chain. Sec. 0 focuses on the dynamics of quantum phase transitions. 
We explain there how fidelity can be used to obtain the probability of finding the system in the ground state after the 
non-equilibrium crossing of the quantum critical point. Finally, we briefly discuss in Sec. IVII two problems associated 
with the central spin coupled to the Ising chain. Their solution also involves fidelity. 


II. QUANTUM ISING MODEL 

The Hamiltonian of the one dimensional quantum Ising model in the transverse field reads 

N 

77(5) = -EK<+i+5<). (8) 

i=l 

where g stands for an external magnetic field acting on the system and N is the number of spins. The periodic 
boundary conditions are assumed, i.e., = erf- While the spin-spin interactions in this model try to align the 

spins in the Fx direction, the magnetic field tries to put the spins along its direction, i.e., +z for g > 0 and —2 
for g < Q. The competition between different spin orientations results in the quantum phase transition between the 
ferromagnetic and paramagnetic phases (Fig. [T]). The system is in the ferromagnetic phase for |g| < 1, it is in the 
paramagnetic phase for |g| > 1, and it is at the quantum critical point when |( 7 | = I, i.e., gc = ±1- The experimental 
possibilities for the studies of this model and other Ising-like models appear in cobalt niobate E3) ‘^old 

ion E3,[il, cold atom [l^, and NM R simulators of spin systems. 

This model is exactly solvable [2ll - [^ . Its role in the studies of strongly correlated quantum systems parallels the 
role of the harmonic oscillator in the single-particle quantum theory: Numerous theoretical concepts in many-body 
quantum physics are being tested in the quantum Ising chain in the transverse field. 

To correctly compute fidelity of the quantum Ising model, one has to properly identify the parity symmetry of 
its ground state [a, [^. We will follow below the discussion from Ref. (^, where the comprehensive study of this 
problem in a finite Ising chain is presented (see also Ref. [2^ for the discussion of the finite-size effects in the XY 
model). 

First, we note that the Hamiltonian ([8]) commutes with the parity operator 

N 

[H,p]=o, 

i=l 

whose eigenvalues are either -fl or —I. Thus, the Hamiltonian (|5]) can be independently diagonalized in the positive 
and negative parity subspaces. 
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Figure 1: Schematic illustration of the quantum phase transition in the quantum Ising model. Adapted from Ref. 


Second, one performs the Jordan-Wigner transformation mapping the spin-1/2 operators to fermionic operators 


< = 1 - 2clc„ < = 


Substituting this into Eq. (|8]), one obtains 


^ (ci + c|) n(^ “ ^ ’ ^ 


j<i 


N 

Hid) = -Yl [/*.*+! + 9(c^cl - cjci) 


2=1 


( 9 ) 


fid = cjcf - Cic] - CiCj + cM, 


where 


CN+1 = —Cl 


( 10 ) 


should be taken in the positive parity subspace and 


Cn+1 = Cl 


( 11 ) 


in the negative parity subspace. Since the Hamiltonian ([S]) is quadratic in the fermionic operators, it can be diago¬ 
nalized in the standard way. One first substitutes 


c, = ^^P( ^/^) exp(iA:j) 

if=±fc 


\/F 


( 12 ) 


into the Hamiltonian ©, where [Hi 
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in the positive parity subspace and 


k = 0, 

fc = 0, 


27r 47r 

)/V’TV’ 

27r 47r 

)/V’TV’ 


, TT —— for even TV, 
TV 

, TT for odd TV, 


, TT for even TV, 

TT 

, TT-- for odd TV, 

TV 


(13) 


(14) 


in the negative parity subspace. The quantization of momenta (1131) and (|14l) follows from the boundary conditions 
m and m, respectively. Then one performs the standard Bogolubov transformation to diagonalize the resulting 
Hamiltonian. 

Third, one needs to determine whether the ground state of the Ising Hamiltonian ([8]) lies in the positive or negative 
parity subspace. This can be done by computing the lowest eigenenergy in each of the subspaces, say £+ and e~ in 
the positive and negative parity subspaces, respectively. We find that e~ — e+ equals 


r,N 


dt 4TVf;;-3/V(l-0(l-'7'l) 


1 - 
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in the ferromagnetic phase, 


sig.(g»)(2|s| - 2) + f di— 


1 — /g'^^ 


in the paramagnetic phase, and 


sign(gf) 


at the critical points. 

A quick look at these expressions shows that the parity of the ground state of the Hamiltonian ([5|) is 


• positive for all magnetic fields when the number of spins N is even, 

• positive for g > 0 and negative for g < 0 in the odd-sized systems. 


These remarks follow from the observation that sign(£“ — £+) = sign((;^). The same result was obtained in a different 
way in Ref. Q for g > 0. 

Besides providing the parity of the ground state, these expressions quantify the gap between the eigenenergies of 
the positive and negative parity ground states of the Hamiltonian (|5]): 

• Ferromagnetic phase (|(;| < 1). The gap is exponentially small in the thermodynamically large systems 

|£“ — £+1 = O(^exp{—N/^{g)) /'/N^ , 

where 

C(5) = l/Nff|| (15) 

is the correlation length of the infinite Ising chain and we define the thermodynamic limit by the relation 

N » e(g). (16) 

In the opposite limit, 

N « C(5), 

we hnd that the gap is polynomial in the system size 


£ — £ 


= 0{1/N). 


• Paramagnetic phase (j^l > 1). The gap is always macroscopic and grows approximately linearly with the distance 
from the critical point. 

• Critical points = 1. The gap is inversely proportional to the size of the system in large systems. 

At the risk of stating the obvious, we mention that the qualitative content of these remarks is discussed in textbooks 
on quantum phase transitions (see e.g. Ref. Q)- 


III. FIDELITY SUSCEPTIBILITY OF QUANTUM ISING MODEL 

Our discussion in this section will be based on Refs. [2^1 and M- The former one provides the exact closed- 
form expression for fidelity susceptibility of even-sized Ising chains and studies its basic properties. The latter one 
generalizes this result to the odd-sized chains and discusses the pitfalls of incorrect identihcation of the parity of the 
ground state in hdelity susceptibility computations. 

The ground state wave-function of the Ising chain can be written as 


\g) =n(‘=°®(Y I 


sin(6»fe) = 


in(y) clclfe) |vac), 
sin(fc) g — cos(fc) 


- 2gcos{k) + 1 


cos(6'fe) = 


\/g^ - 2gcos{k) -f l’ 


(17) 
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Figure 2: Panel (a): black solid line shows (1191) . while the red dashed line shows x~ II20II . Even = 40 is used. Panel 
(b): The X^ jxT ratio for data from the (a) panel. Panel (c): the same as panel (b), but as a function of the ratio between the 
system size and the correlation length Nl^{g) = N\ ln|^| |. 


where the proper momenta k have to be used: (TT^ if the ground state of the Hamiltonian (|5]) lies in the positive 
parity subspace and (IT^ if it belongs to the negative parity subspace (see Ref. (U for the discussion of the k = 0 ,tt 
modes that do not contribute to the following calculations). The |vac) state is annihilated by all Ck operators. 

Using this wave-function one easily finds with Eq. ([S]) that 


xig) = 



1 sin^(A:) 

^ k -‘^9 cos{k) + if 


(18) 


This summation has been analytically performed in Refs. and (^ . 

Before providing the definite closed-form expression for fidelity susceptibility, we emphasize the need for the correct 
identification of the parity of the ground state 0 . If we use the momenta (HU to compute the series m, we will 
obtain 


X^(ff) = 


iV2 




N 


9^ -g"^ 


16g'^[gN^.l f 16g2(^Ar X)(5f2 _ !)■ 


If we, however, sum the series (fTSl) over momenta da, we will get 


X (5) = - 


iV2 




N 


5^ + 5^ 


IQg^ig^-lf mgfg^ -W-1)' 


(19) 


( 20 ) 


These two expressions are plotted in Fig. [5^. They differ negligibly away from the critical point, i.e., for magnetic 
fields g such that N ^(g). On the other hand, near the critical point, i.e., for magnetic fields g such that the system 
size N is smaller than a few correlation lengths ^( 5 ), they substantially differ (Fig. [^). This difference was numerically 
studied in Ref. where it caused problems in the Quantum Monte Carlo studies of fidelity susceptibility at non-zero 
temperature. Our expressions analytically quantify this difference in the zero temperature limit. For example. 


• right at the critical point gc = 1 the ratio of the two expressions equals 


X“ 


= 3 


N 




IV- 2’ 


and so it approaches 3 as N increases (Figs. [5 Jd and [5}:), 


















6 



magnetic field g 


Figure 3: Black solid line shows x+ (HU. while the red dashed line shows x ■ (unj. The plot is done for odd A'’ = 41. 


• in odd-sized systems neither of the two expressions, (HU and (EOl), is symmetric with respect to the g O —g 
transformation (Fig. [3]). 

Using the proper identification of the parity of the ground state discussed in Sec. m one quickly finds that (i) 
x{9) = X^{g) for 8.11 system sizes and g > 0; (ii) xig) = X^[9) for even-sized systems and g < 0; and (hi) x(g) = x~(9) 
for odd-sized systems and g < 0. Correlating the observations (i) - (iii) with equations (fT31) and (I3D1) . we get 


N \g\^ — g^ 

= W{\gr + lf ^ l6gH\gr + l){g^-iy ^ 

This exact closed-form expression works for any system size N > 2 and any magnetic field g. As expected, it is 
symmetric with respect to the g o —g transformation. Its properties can be summarized as follows 
First, it undergoes a simple transformation law under the Kramers-Wannier duality mapping 


Indeed, one easily verifies that 


1 

9 ^-. 

9 


9^X{9) 



( 22 ) 


(23) 


Second, one finds that away from the critical point in the paramagnetic phase 

- 1 ) + 

where the remainder R, which can be easily computed from Eq. m, can be shown (2^ to be negligible when 
inequality (1161) holds. It should be stressed that the main part of this result can be obtained by replacing the sum by 
the integral (see e.g. Refs. i, 0 , and 0) 



(25) 


in Eq. (1181) . While such a replacement works well away from the critical point, it leads to a completely wrong 
result near the critical point. The formal approach to the replacement of the sum by the integral is provided by the 
Euler-Maclaurin summation formula . For example, adopting this formula to the even N positive parity case, we 
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get 


where 



f{k) = 


1 

4(7 


sin^(fc) 

2gcos{k) + 1)^ ’ 


/(-)(fc) 


dfc'" 


/(fc), 


(26) 


an integer s > 1 is a parameter that can be chosen at will, are Bernoulli numbers, and (j)s are the Bernoulli 
polynomials (we use the convention for Bg and (pg from Ref. (29|). Therefore, by performing (|25l) one skips the 
second, third and the fourth line of the right hand side of Eq. (HH), and replaces the integration ranges ir/N and 
TT — tt/N in the first line by 0 and tt, respectively. This procedure is bound to fail near the critical point. A quick 
inspection of expression (l26l) shows that the exact computation of the sum (ITS)) is easier than the formal study of the 
Euler-Maclaurin expression (1261) ! 

Before moving on, we mention that the application of the duality mapping (1^^ to Eq. ((Ml) provides us with the 
following expression for fidelity susceptibility in the ferromagnetic phase 


xig) = 


N 


Ra/\9\) 


16(1-52) g 4 : 

where again the remainder is negligible away from the critical point, i.e., when the condition (11611 is satisfied. 
Third, near the critical point one can expand the exact expression in a series to get 


(27) 


xig) = 


NjN-l) / 

3252 


+ 1 iN\n\g\)^ 
~ 6 


+ ^(l5l) 


where this time the reminder r(| 5 |) is negligible for N <C C(5)- In particular, at the critical points one finds 


xigc) 


N{N - 1) 
32 


(28) 


Eourth, noting that d = v = 1 and = if in our quantum Ising model (^ . one easily verifies the scaling relations 
® and ([T|) by comparing them to the equations (IMl) . (IMI) . and (1^^ . 

Eifth, one can show from the exact solution (ED), that the maximum of fidelity susceptibility is located in the 
ferromagnetic phase at the distance 


— - — +0(N-'^) 

1V2 ATS ^ I ) 

from the critical point [ 2 ^. The fact that the maximum is shifted away from the critical point due to the finite-size 
effects is not surprising: The shift proportional to 1 is expected for free fermionic systems [s^ • The shift of the 
maximum into the ferromagnetic phase is the consequence of the duality symmetry (I23p , which the reader can easily 
verify. 


IV. FIDELITY IN THERMODYNAMICALLY LARGE ISING CHAIN 

Using wave-function (HD, one easily finds that 

i^(5,^)=n‘ 


cos 


9+{k)-9_{k) 


tan5±(fc) = 


sin(fc) 


g±6 — cos(/c) ’ 


( 29 ) 
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Figure 4: Illustration of the transition to the thermodynamic limit. In both panels the black solid line shows = 1,5)- Upper 
panel: the field shift S — lO”"*. Lower panel: the system size N is set to 10®. Red dots come from F = 1 — N(N—1)S^/16, i.e., the 
approximation of fidelity m through fidelity susceptibility (1281) . The black crosses come from the leading-order thermodynamic 
approximation (EU, which for g = 1, i.e., c = 0, reads InF = —N\5\/A:. 


We are now interested in studying this expression in the thermodynamic limit near the critical point gc = 1. In fact, 
from now on we will focus on g > 0 and assume even N for convenience. This section is based on our Refs, [l^ and 

[H. 

The transition to the thermodynamic limit can be induced by either the change of the system size N or the change 
of the field shift 5 (Fig. [4]). It happens in our problem when 

N > min[(^(5 -b 5),^{g - (5)]. 

Substituting ± d) dT51) . we find that near the critical point this condition implies that the system reaches the 
thermodynamic limit when 


N\S\ » 1, (30) 

which is numerically studied in Ref. [l3|. It explains why the change of either N or 5 can drive the system into the 
thermodynamic limit. 

It is currently unknown how to find a closed-form expression for the product ([29]). Thus, we have to restore to 
approximations. We assume that oo at the fixed S, take the logarithm of Eq. (|29l) . and then replace the sum 

over momenta with the integral (I25p . 

Introducing the relative distance from the critical point 


c = 


g- 1 

l<5| 


we obtain the following approximate result |13l . I31l 


lnF(g,5) 


N 




( 31 ) 
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where 


1 ^ aaO (Id - l)lmE(q) 
fcl lcF/«c,) (Icl-ll-- 


Cl = -4- 


/TT 47r 

A K{ci) _ (|c| - l)Im£;(c 2 ) 

27r 47r 

(|c| + l)^ 

C2 = 


(32) 


■(Id -1)2’ (|C|-1)2' 

The complete elliptic integrals of the first and second kind are defined as 

V2 rfd p. . 

^=, E[x) = 

/o 


K{x) = / — . , E{x) = / dM 1 — X SIT? {(j)), 

yl —a;sin^(d) 


respectively. Several remarks are in order now. 

First, Eq. (1511) provides the leading order term of the TV —>■ oo expansion of the exact expression for fidelity per 
lattice site, i.e., — Inf (g , Sy N. Its accuracy is illustrated in Fig. [SI The next order correction to fidelity per lattice 
site is discussed in Ref. [^, where we refer the reader for the details. It comes from the difference between the sum 
and the integral: 


tt—tt/N 

In cos 

k—ir/N 


0+{k)-0-ik) 


N 

2tt 


dk In cos 


0+{k)-0_{k) 


In particular, the next order approximation to hTF[g,5)/N at the quantum critical point = 1 is 


lnF(l,(5) 

N 


4 


ln(2) 
2N ’ 


(33) 


(34) 


where the ln(2) /2 prefactor of the subleading term comes from the evaluation of Eq. (ISSl) in the N ^ oo limit |31j . 
We will show in Sec. j^that the subleading term can be accurately extracted out from the probability of finding the 
system in the ground state after the non-equilibrium quantum phase transition. 

Second, it should be stressed that Eq. (1511) is derived for |^| <C 1 and Iff — 1| = |c5| <C 1. This is sufficient for 
describing the sharp drop of fidelity in the vicinity of the critical point. It should be noticed that fidelity in the 
thermodynamic limit can change by several orders of magnitude around the critical point when the magnetic field is 
changed by a few S’s. Note that this observation holds for 1/TV <C |^| 1, so fidelity is indeed a sensitive probe of 

quantum criticality! 

Third, Eq. (|31l) analytically shows how the Anderson orthogonality catastrophe happens in the quantum Ising 
model. It proves that the overlap between the Ising model ground states in the quantum critical region is exponentially 
small in the system size in the thermodynamic limit. 

Fourth, the expression for A{c) can be simplified for |c| 1. Indeed, one finds that A{c) approaches l/16|c| in this 

limit (Fig. O. This leads to the following approximation 


F{g,6) ~ exp - 


N5^ 




(35) 


If additionally —N6^/16\g — 1| <C 1, then this expression yields the same result as Eq. ©. Finally, we note that the 
solutions (l?T|) and (1351) agree with the scaling results ([5]) and ([7]), respectively. 


V. DYNAMICS OF QUANTUM PHASE TRANSITIONS 


We will discuss now what insights into the dynamics of the quantum quench are provided by fidelity. Although 
the discussion below is based on Sec. V of Ref. [3l|, the numerical results presented here have not been published 
before. We focus here on the dynamics of the quantum Ising chain, but the discussion can be easily generalized to 
other “typical” quantum critical systems (see the end of this section). 

To start, we consider an instantaneous quench, i.e., a sudden change of the parameter g of the Hamiltonian (|8]) 
from gi to g 2 - If the system was initially prepared in the ground state |gi), then the probability of finding it in the 
ground state after the quench equals 


\{9i\92)\=F^ig,6), g = 


9i +92 


6 = 


92 - 9i 


2 


2 
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c 


Figure 5: The scaling function A{c) (I32II . Black line shows —A{c) while the red crosses show ln_F(l + c|(5| ,5) /Af|(5|, where 
N = 10® and <5 = 7r/1000. Note that the thermodynamic limit condition is well-reached for such a system as = IOOtt ^ 1. 
Adapted from Ref. 0. 
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Figure 6: Illustration of the adiabatic-impulse approximation in the quantum Ising chain. 


where F{g, 6) is given by Eq. (1^^ and the approximations discussed in Secs. IIIIl and lTVl can be used in the appropriate 
parameter ranges. 

A much more interesting quench dynamics shows when the magnetic field is continuously changed in time 

9{t) = —(36) 
TQ 

where t : —oo —>■ 0 and the quench time tq controls the speed of the driving. The dynamics of the Ising chain due to 
such time variation of the magnetic field is exactly solved in Ref. [s^l (see Refs, [sj and [sj for the recent reviews of 
the dynamics of quantum phase transitions). We will now link fidelity to the probability of finding the system in the 
ground state after such a quench. To do so, we invoke the quantum version of the adiabatic-impulse approximation 
discussed in Ref. (see also Ref. M)- This approximation splits the evolution across the critical point into three 
parts depicted in Fig. [HI 

• The first adiabatic regime takes place when the system is far away from the critical point, say 

g(t) >l + g, 

where g will be discussed below. It is assumed in this regime that the system’s wave-function |'0(t)) ~ l5(0)i 

where \g{t)) denotes the instantaneous ground state of H{g{t)) ([5]) and is the global phase. Such dynamics 
is possible when the gap in the excitation spectrum is large enough. 
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Figure 7: Panel (a): the probability of finding the N = 150 Ising chain in the instantaneous ground state during a quench. 
Panel (b): the probability pcs of finding the Ising chain in the ground state by the end of time evolution. The quench protocol 
in both panels is provided by Eq. dMli with TQ = 50. The evolution starts from the paramagnetic phase ground state at the 
magnetic field g(t = —5tq) = 5 and it ends in the ferromagnetic phase at the magnetic field g{t = 0) = 0. The solid line on 
both panels is a numerical calculation. The red dashed-dotted line in the lower panel shows 1 — exp(—27r®rQ/A^^). It illustrates 
what happens when the finite-size effects play a role (3^ . 


• The impulse regime happens around the critical point, 

\g{t) - 1| < 

and the wave-function, apart from the overall phase factor, is assumed to be constant: \tp{t)) ~ -I- g). 

This happens when the gap in the excitation spectrum becomes too small to support adiabatic dynamics at the 
given quench time. 

• The second adiabatic regime happens when 


9{i) < 1 - 

where it is assumed that no additional excitations are created by the quench, i.e., \{g{t)\ilj{t))\ = const. 

These three stages of evolution are easily seen on Fig. [7^, where the squared overlap of the wave function |V'(t)) onto 
the instantaneous ground state \g{t)) of the Hamiltonian (|5]) is plotted for a typical time evolution. We mention in 
passing that we have done the numerics for this section in the same way as in our recent Ref. (3^ discussing the 
counterdiabatic dynamics of the Ising chain j^. 

The size of the impulse regime is governed by the gap in the excitation spectrum and the rate of driving. It is 
estimated by solving the following equation (see e.g. Refs, [s^, [s^, and for different versions of this simple 
estimation and the physical reasons supporting it) 

1 _ gap 

gap(l±5) ^ ^gap 

Assuming that the system is infinite, the relevant energy gap in the quantum Ising model is given here by 2\g — 1| 
for g > 0 that we consider in this section. This gives us 


1 


9 


(37) 
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Figure 8: Logarithm of the probability of finding the Ising chain in the ground state after the quench; see the discussion around 
Eq. (I42II for the details. 


It should be stressed that there is no sharp splitting between the adiabatic and impulse dynamics. Instead, there is 
a crossover and Eq. (EZD provides the scaling of its location with the quench time tq . Such estimation is meaningful 
for slow quenches only, rg ^ 1, which can be numerically checked. 

The quench time, however, cannot be too large in a finite Ising chain if we want to see the impulse regime. Indeed, 
given the fact that the gap at the critical point scales as 1/iV, one can always drive the system so slow that it will 
adiabatically cross the critical point (its excitation will be then exponentially small in the quench time [s^, which 
we illustrate in Fig. IZb). Note that the gap in this section refers to the difference between the eigenenergy of the 
ground state and the lowest excited state that can be populated during the evolution. It shall not be confused with 
the gap between the positive and negative parity subspaces studied in Sec. im the parity is conserved during the time 
evolution. 

We expect that the finite-size effects are negligible during a non-equilibrium quench, when the spin-spin correlations 
in the system’s wave-function \ijj(t)) decay on a length scale much smaller than the system size. Within the adiabatic- 
impulse approximation this length scale is upper bounded by the equilibrium correlation length (|15l) taken at 1 -\- g. 
Thus, it is reasonable to expect that the finite-size effects are negligible when 

(38) 

which we assume below. This educated guess is rigorously discussed in Refs. and (s^- 

Using the adiabatic-impulse approximation, we estimate that the probability pcs of finding the system in the 
instantaneous ground state away from the critical point, i.e., when |g — 1| ^ g, is 

PGS = \{l-9\l + 9)?=F^{l,5), (39) 

5 = 9^^. (40) 

Combining Eqs. (1551) and (HOI) , we see that the condition (1501) is satisfied. Therefore, the thermodynamic expression 
(154)) for fidelity should be used in Eq. (1591) . which leads us to 


PGS = 2 exp 



const 


( 41 ) 


where the prefactor of 2 comes from the subleading term in Eq. (1551) . 
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This simple prediction can be compared to numerics. To do so, we employ the quench protocol (1361) . The evolution 
starts from the paramagnetic phase ground state at the magnetic field g(t = —5tq) = 5 and it ends in the ferromagnetic 
phase at g(t = 0) = 0. To simplify the comparison between the theory and numerics, we take the logarithm of Eq. 
(liTj) getting 


In(pGs) = ln(2) - N X 
ln(2) = 0.693147-•• . 


(42) 


First, we choose the quench time tq = 50 and change the system size N in the range [100,1000]. The standard 
linear ht to the numerics from Fig. [5^ yields: 

lii(pGs) = 0.693146(1) - 0.020800730(1)7V. 

The linear dependance of In(pGs) on the system size TV is nicely confirmed by the fit, and we note that the fitted 
intercept surprisingly well agrees with the theoretical one. 

Second, the system size is set to TV = 150 and the quench time tq is varied in the range [50,150]. The linear fit to 
the numerics from Fig. [SId gives 


. X . X 22.065(1) 

In(pGs) = 0.6940(1)-——. 

The intercept again agrees surprisingly well with the theoretical prediction and the linearity of In(pGs) in th® inverse 
square root of the quench time is nicely confirmed. 

The fitted lines are not depicted on Figs. [8^ andlSh because they are practically indistinguishable from the numerics. 
It is worth to stress that for the system sizes considered in these numerical studies the inclusion of the subleading 
term to fidelity per lattice site dMl) is crucial for getting a good agreement between the theoretical prediction for pcs 
and numerics. 

Finally, we mention that a more general result can be obtained if we assume that the gap closes at the critical point 
in the infinite s yst em as \g — gc\^‘^, where z and v are the universal critical exponents. Repeating the steps outlined 
above one gets 


PCS 



const ^ 
dv/{l+zu) I ’ 
'^Q / 


(43) 


where d stands for the dimensionality of the system. The Ising result is recovered for z = ly = d = 1. 

The results (HTl) and (H51) have simple interpretation in the context of symmetry breaking phase transitions, which 
in the quantum Ising model happen when the system is driven from the paramagnetic to the ferromagnetic phase. 
As predicted by the Kibble-Zurek theory [4l|, such transitions lead to the creation of topological defects whose 
density scales with the quench time as Thus, the probability of finding the system after the quench in 

its instantaneous ground state is exponential in the number of topological defects created during the non-adiabatic 
crossing of the quantum critical point. 


VI. FIDELITY IN CENTRAL SPIN SYSTEMS 

It turns out that knowledge of fidelity helps in the studies of the central spin-1/2 coupled to the Ising chain (4^1^. 
The Hamiltonian of such a system reads 


N 

1=1 

where H(g) is given by Eq. ([5]), is the Pauli matrix of the central spin, and d is the central spin-environment 
coupling. Such a model is known as the central spin model because of the uniform coupling of the central spin to 
its “environment”. Without going into details, we will mention two problems, where the expressions for fidelity from 
Secs. Ell and 113 can be used. 
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Critical dynamics of decoherence. Suppose that the magnetic field is quenched by the protocol (1361) . The system is 
initially in the product state 


{c+\t)+c-\l))(^\g{t =-oo)), (44) 

where the central spin is in an arbitrary superposition of its up/down states and the Ising chain is in the ground 
state \g) of H{g). The idea now is to study the decoherence of the central spin due to the presence of the many- body 
Ising environment undergoing a non-equilibrium quench across the critical point (s). As is explained in Ref. |42l |. 
the knowledge of fidelity F{g,6) is crucial for the explanation of the decoherence rate of the central spin in this 
non-equilibrium system. 

Magnetic Schrodinger’s cats. We consider now the same system as the one described above and assume that the 
Ising chain is adiabatically driven towards the critical point, where the time variation of the magnetic field (I36|) stops. 
Assuming that the evolution starts from the initial state (IT41) and that tq ^ N'^, the final state of the system will be 
approximately given by 

14.) (^\g-S) + c+e*'^+ It) O Iff + S), 

where (f>± stand for the phases picked by the system during the evolution [4^. Assuming that g — S<gc<g + S, 
Iff — 5) and Iff + 5) are the ferromagnetic and paramagnetic ground states of the Ising Hamiltonian. This state is the 
magnetic Schrodinger’s cat state. The dead/alive cat states are the ferromagnetic/paramagnetic ground states. The 
state of the radioactive atom that has/has not decayed is played by the up/down states of the central spin. It is not, 
however, the ordinary Schrodinger’s cat because the two cat states are not orthogonal, 

\{g-S\g + S)\= Fig, 6)^0. 

It would be now certainly interesting to understand how the rate of decoherence of such a superposition state depends 
on the similarity between the two cat states, i.e., their fidelity. It is perha ps w orth to stress that the experimental 
studies of the decoherence of “simpler” superpositions have been recognized [31^ by the Nobel committee in 2012, 
so it would be certainly very exciting to experimentally quanti fy th e decoherence rate of the magnetic Schrodinger’s 
cats! For other applications of fidelity in this context see Ref. [^. For similar ideas in the cold ion and cold atom 
systems see Refs. [4^ and respectively. 


VII. SUMMARY 

We finish these lecture notes by quoting P. W. Anderson (l^ “ While wave functions and overlap integrals are often 
of little consequence in many-body systems, this one is at least related to the response to a sudden application of the 
potential, and indicates that response involves only the emission of low-energy excitations into the continuum, as well 
as that the truly adiabatic application of such a potential to such a system is impossible. We hope that these lecture 
notes provide further evidence that the wave functions and overlap integrals are indeed very useful in the studies of 
many-body quantum systems. 
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